Reverse epitaxy of Ge: ordered and facetted surface patterns 
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Normal incidence ion irradiation at elevated temperatures, when amorphization is prevented, 
induces novel nanoscale patterns of crystalline structures on elemental semiconductors by a reverse 
epitaxial growth mechanism: on Ge surfaces irradiation at temperatures above the recrystallization 
temperature of 250° C leads to self-organized patterns of inverse pyramids. Checkerboard patterns 
with fourfold symmetry evolve on the Ge (100) surface, whereas on the Ge (111) surface, isotropic 
patterns with a sixfold symmetry emerge. After high fluence irradiations these patterns exhibit well 
developed facets. A deterministic nonlinear continuum equation accounting for the effective surface 
currents due to an Ehrlich-Schwoebel barrier for diffusing vacancies reproduces remarkably well our 
experimental observations. 



PACS numbers: 68.35. Ct, 79.20.Rf, 81.16.Rf, 81.65.Cf 

Self-organized pattern formation in systems far from 
equilibrium is a fundamentally interesting phenomenon 
governed by the interplay of kinetic and diffusive mech- 
anisms. In these systems the dynamics of pattern for- 
mation obeys universal laws with general scaling rela- 
tions pQ. In addition, surface patterns with nanoscale 
dimensions are of technological interest for applications 
in sublithographic surface templating and for quantum 
dot device fabrication [2 j. They may be generated on 
surfaces by homoepitaxy [3] , heteroepitaxy [4 or by ener- 
getic ion irradiation [5H7]. However, semiconductor sur- 
faces become amorphous during conventional ion irradi- 
ation at room temperature. At these conditions periodic 
ripple patterns oriented perpendicular or parallel to the 
ion beam direction and isotropic, hexagonally ordered, 
dot or hole patterns have been observed [7] and are inde- 
pendent of the crystal structure [8]. The origin of these 
patterns is attributed to an interplay of a surface insta- 
bility due to sputtering and mass redistribution together 
with surface relaxation mechanisms [9j [10} . Below inci- 
dence angles of 50° smoothing is expected to dominate 
on amorphized elemental materials [11]. On the other 
hand, metals and metal oxide surfaces remain crystalline 
during ion irradiation at room temperature [6] . They ex- 
hibit a much higher complexity of pattern formation due 
to additional instabilities resulting from anisotropics in 
surface diffusion and due to biased diffusion across step 
edges [12]. Strong similarities with homoepitaxy have 
been identified: layer by layer erosion has been observed 
on metals [13] as well as on semiconductor surfaces [T4] . 
Furthermore, similar to mound formation in epitaxy, pit 
formation has been observed experimentally on ion irra- 
diated metal surfaces [15j [16] . Although the formation 
of pits has also been seen in low fluence irradiations of 
semiconductors [T7HT9] , dense and ordered patterns of 
facetted nanostructures, as found in homoepitaxy and 
heteroepitaxy, have not been observed until now on ion 
irradiated semiconductor surfaces. 



In this letter we present the formation of regular pat- 
terns of crystalline structures induced by normal inci- 
dence ion irradiation of an elemental semiconductor, Ge, 
at elevated temperatures. At temperatures above the re- 
crystallization temperature of 250° C, ion induced bulk 
defects are dynamically annealed and the surface re- 
mains crystalline. By only increasing the irradiation tem- 
perature, instead of inducing surface smoothing, novel 
checkerboard patterns with crystalline facets appear. 
These patterns exhibit the symmetry of the crystalline 
surface, i.e. on the (100) surface the structures are ori- 
ented along the (100) direction whereas on the Ge (111) 
surface pit patterns with a six fold symmetry develop. 
They strongly resemble mound patterns in homoepitax- 
ial and heteroepitaxial growth [3 , but are reversed. The 
mechanism can thus be interpreted as reverse epitaxy. 
In analogy with the Villain instability in homoepitaxy 
resulting from the Ehrlich-Schwoebel (ES) barrier for 
an adatom descending a monoatomic step [20j [21] we 
conclude that the formation of ion induced crystalline 
patterns results from the existence of an ES barrier for 
the ascending of surface vacancies created by sputtering. 
Based on the proposed atomistic mechanisms we derived 
a continuum equation for reverse epitaxy which describes 
remarkably well the experimentally observed surface evo- 
lution. 

Ge (100) and Ge (111) samples cut from epi-ready Ge 
wafers (root mean square roughness w rms = 0.7 nm) were 
irradiated by a broad 1 keV Ar + ion beam at normal 
incidence without any pre-treatment. The irradiations 
were performed in a high vacuum chamber (base pres- 
sure 10 -8 mbar) equipped with a Kaufman ion source 
with a single graphite grid extraction with 5 cm diam- 
eter. During irradiation the chamber is backfilled with 
Ar at 10 -4 mbar. In order to avoid metallic contamina- 
tions from the sample holder the samples with a typical 
size of 5 x 5 mm 2 were glued to a Si plate with size of 
10 x 15 mm 2 . The samples were heated by a boron ni- 
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FIG. 1. AFM height images (4x4 \im 2 ) of Ge (100) sur- 
faces after ion irradiation with an ion fluence of 3 x 10 18 cm -2 
at temperatures of 230°C, 260°C, 350°C, and 430°C, respec- 
tively. The (100) crystal directions are marked by arrows in 
(b) and (c). The height color bar on the right side represents 
8 nm for (a) and (d), and 38 nm for (b), (c), respectively. 



tride heater from the backside of the Si plate. The sur- 
face temperature was monitored during ion irradiation 
by a pyrometer operated in the wavelength range of 2- 
2.8 /im that was calibrated externally by a thermocouple. 
The surface topography was analyzed after irradiations 
ex-situ by a multimode atomic force microscope (AFM) 
from Veeco in tapping mode. 

In Fig. [I] AFM images of Ge (100) surfaces irradi- 
ated with ion fluences of 3 x 10 18 cm -2 and different 
surface temperatures ranging from 230° C to 430° C are 
shown. At temperatures below 250° C the Ge surface 
remains smooth exhibiting no patterns after irradiation 
(Fig. [TJa)). At these conditions the Ge surface is amor- 
phized by ion irradiation and smoothing by surface diffu- 
sion and mass redistribution dominates [11]. The initial 
roughness, w rrnsi of the virgin Ge surface of 0.7 nm is 
reduced to 0.18 nm. At temperatures higher than 250° C 
checkerboard patterns appear after irradiation indicating 
that an additional ion induced instability appears. The 
structures of the pattern have a rectangular shape and an 
average size of ~ 150 nm with an orientation in the (100) 
direction (Fig.Qb)). The size of the structures increases 
with temperature. A checkerboard pattern with an av- 
erage size of ~ 260 nm develops under ion irradiation at 
350° C (Fig. Qc)). Again, the pattern is oriented along 
the (100) direction. At 430° C the structure size increases 
further while the symmetry is changed to an isotropic 
pattern of pits with diameter of 300 — 1000 nm. The 
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FIG. 2. TEM images of Ge surfaces in cross section irradi- 
ated at temperatures of (a) 230°C, (b) 260°C, and (c) 350°C, 
respectively. In (a) a 2.6 nm thick amorphous layer can be 
identified on top of the crystalline bulk, whereas in (b) and 
(c) no amorphous layer is visible. 



order is much lower and the structure size distribution 
much broader than for the checkerboard patterns. Fi- 
nally, at temperatures above 500° C the surface is again 
smoothed by ion irradiation, similar to irradiations at 
230° C (not shown). 

We investigated the microstructure of the patterns 
in cross section with transmission electron microscopy 
(TEM). In Fig. [2] TEM images of Ge (100) surfaces are 
shown, which were irradiated at temperatures of 230° C, 
260° C, and 350° C, respectively (same irradiation con- 
ditions as in Fig.[l{a)-(c)). At 230°C the surface is flat 
and a 2.6 nm thick amorphous Ge layer is visible, in good 
agreement with the calculated range of 3 nm for 1 keV 
Ar + ions in Ge [22]. At 260° C and higher no amorphous 
layer is visible. In the high resolution TEM images of the 
Ge surfaces irradiated at 260° C and 350° C flat facets are 
visible. The angle between the facets and the (100) plane 
is determined to 8°-10°. 

The transition from smoothing to roughening by ion 
irradiation between 230° C and 260° C can be clearly at- 
tributed to the temperature at which amorphization by 
ion irradiation is prevented. Ion induced vacancies and 
interstitials in the bulk are dynamically annealed above 
this temperature and only adatoms and surface vacancies 
remain as defects. Hence, surface patterns result from 
the kinetics of adatoms and vacancies created by the in- 
cident ion beam, however, vacancy kinetics is expected 
to dominate, because the number of created vacancies is 
higher than the number of adatoms by sputtering. In 
analogy to the ES barrier for adatoms, i.e. the diffusion 
barrier for an adatom to descent a step edge, a barrier 
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FIG. 3. AFM images with size of 500 x 500 nm of (a) Ge (100) 
surface and (b) Ge (111) surface irradiated at temperatures 
of 280°C. 



exists for a vacancy to ascend into the next higher ter- 
race. Thus, vacancies are trapped on lower terraces lead- 
ing to 3D reverse growth of surface structures, i.e. the 
formation of pyramidal pits (inverse mounds). At higher 
temperature the mobility of vacancies increases and the 
nucleation density decreases leading to pits with larger 
separation and size. If the thermal energy of vacancies is 
high enough to overcome the ES barrier, the surface will 
remain smooth which is in agreement with our observa- 
tion that above 430° C the surface remains flat. 

According to the terraces-step-kink (TSK) model the 
evolution of crystalline surfaces is described by atomistic 
processes on terraces, steps, and kinks [23]. Diffusion 
and attachment of adatoms and vacancies at steps and 
kinks is expected to be anisotropic on crystalline sur- 
faces. Therefore, the pyramidal pits exhibit facets ori- 
ented along the crystalline directions. On Ge (100) the 
edges of the pits are oriented along the (100) crystal di- 
rections, whereas on Ge(lll) patterns with a six fold 
symmetry appear (see Fig. [3|. The emergence of these 
edges is attributed to an additional barrier at kink sites, 
similar to the barrier at step edges [24]. On the (100) sur- 
face the kink ES barrier is responsible for the repulsion 
of vacancies on the fast diffusing (110) step edges. This 
mechanism has also been proposed for the formation of 
mounds oriented in the (100) directions in homoepitaxy 
of Ge (100) [25] . The energy barrier for the diffusion 
around corners is expected to be smaller than for cross- 
ing step edges. Therefore, the kink ES barrier vanishes 
already at temperatures where the terrace ES barrier is 
still active and isotropic pit patterns should appear. This 
is indeed observed for irradiations at 430° C (Fig. [ljd)) 
where the square symmetry of the pattern disappears and 
dense round pits are formed. 

In order to further elucidate the formation mechanism 
we investigated the roughening and coarsening behav- 
ior of the checkerboard pattern on Ge (100) at 350° C. 
Fig. |4ja) shows the evolution of Ge(100) surface as well 
as the corresponding 2D FFT and 2D angle distribution 
as a function of ion fluence. The FFT reveals the sym- 
metry and the order of the pattern, whereas from the 
2D angle distribution the formation of dominant crystal 
facets can be deduced. After irradiations with a fluence 




FIG. 4. (a) AFM images of Ge (100) surfaces irradiated 
at 350°C and fluences of 1 x 10 17 cm" 2 , 1 x 10 18 cm" 2 , and 
1 x 10 19 cm -2 , respectively. The arrows indicate (100) direc- 
tions in the surface, (b) Snapshots of the numerical integra- 
tion of the continuum equation for 400 s, 2000 s, and 10000 s, 
respectively. The parameters used for the numerical integra- 
tion are e = 1, « = 4, a = — 1, and 5 = 25 on an integration 
grid of 400 x 400 points with Ax = Ay = 1, and At = 0.01s. 
Above the height images the 2D FFT (left) and the 2D angle 
distribution (right) are shown. 



of 1 x 10 1T cm -2 a pit pattern with no clear orientation is 
visible. The slope distribution is isotropic with a maxi- 
mum at 0°. At higher fluence of 1 x 10 18 cm -2 the pattern 
already exhibits an orientation along the (100) crystal di- 
rection. The angle distribution reveals a fourfold symme- 
try with maxima around 9.5°. Finally, at 1 x 10 19 cm -2 
the angle distribution has four distinct peaks around 11°. 
The peaks in the angle distribution are a clear signature 
for facet formation on the patterned Ge (100) surface. 
The azimuthal and polar orientation of the facets can 
be identified with {106} (9.6°) and {105} (10.5°) crys- 
tal planes. The angle is determined by the formation of 
thermodynamic metastable low index planes, which are 
often observed in Ge homoepitaxy and heteroepitaxy on 
Si [2j [4]. The energy minimization of the (105) facet 
results from the effective dimerization of the dangling 
bonds on the two- atom wide (100) terraces separated by 
(100) steps. 

For the description of the temporal evolution of the 
surface during ion irradiation a continuum equation is 
derived considering sputtering and mass distribution by 
the incident ion beam [TJ [9j [10] [26] and surface diffu- 
sion of ion induced vacancies on crystalline surfaces [27] . 
Redistribution and diffusion processes are mass conserv- 
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ing and can therefore be described by surface currents. 
The temporal evolution of the surface height, h(x,y,t), 
can be described by the following deterministic partial 
differential equation: 



dh 



-V + vV 2 h - Vjion - Vjdiff, 



(1) 



where v Q is the constant erosion rate of the flat Ge sur- 
face, vV 2 h the so-called negative surface tension with 
v < 0, which takes into account the curvature dependent 
sputter rate [9 , j{ on the surface current resulting from 
the ballistic mass redistribution induced by momentum 
transfer from ions to surface atoms [10 , 28J, and jdiff the 
surface current due to diffusion [20]. Mass redistribution 
by ion impact is also proportional to the surface curva- 
ture, V 2 /i, with a positive coefficient and dominates the 
destabilizing sputtering term, vV 2 h, at normal and small 
angles of incidence [11] . 

On amorphous surfaces, diffusion is isotropic and is 
described by the Herring-Mullins (HM) surface diffusion 
current [29]. On crystalline surfaces, the diffusive cur- 
rent has to include atomistic surface currents on terraces, 
across terrace steps, along steps, and across kinks [30] . 
On the flat terraces of Ge(100) and Ge(lll) an isotropic 
surface diffusion has to be assumed as well described by 
the HM surface diffusion, j'hm, whereas diffusion across 
steps is biased by the ES barrier for ascending vacan- 
cies, resulting in a net uphill mass current, j'es- The 
symmetry of the crystalline Ge(100) surface is taken into 
account by an anisotropic current vector, Jes'(^), which 
is a function of the surface slopes, rri( x ^ — d( x ^ y \h, in x 
and y direction [31]: 



jdiff = ^V(V 2 /i) + aV(yh) 2 + e 



m x (l 
m y (l 



5m 2 x ) 
5m 2 y ) 



(2) 



In Eqn. 2 the first term on the rhs is the isotropic 
HM surface diffusion, j'hm? the second term is called 
"conserved Kadar-Parisi-Zhang (KPZ)" term and has 
been introduced as a nonlinear current corresponding 
to the "non-conserved" nonlinearity in the KPZ equa- 
tion (A/2(V/i) 2 ) [20]. This nonlinear current is known to 
break the up/down symmetry which is prominently seen 
in the high temperature round pit patterns (Fig. [TJd)). 
However, also the checkerboard patterns in Fig. [TJb)-(c) 
are not symmetric under the transformation h —> —h. 
For positive values of a the surface evolves to patterns 
with mound structures, whereas for negative values of 
a patterns of pit structures appear. Finally, the third 
term describes the anisotropic ES surface current, jes? 
with e the proportionality factor between the ES cur- 
rent and the surface slope [32]. For positive e the ES 
current is pointing uphill inducing a surface instability 
[20] . The anisotropy of the surface current has its ori- 
gin in the anisotropy of the ES barrier itself as well as 
in additional currents along step edges due to step edge 



diffusion [33 , 34 . It is well known that this kind of sur- 
face current leads to the formation of facets at angles 
for which the current becomes zero [27j [30]. Close to 
these points the surface current is negative (positive) for 
smaller (larger) angles, leading to an increase (decrease) 
of the slope. The parameter S determines the angle of 
the facets: = ± arctan^l/^). 

The pattern formation at normal incidence is thus 
dominated by the diffusive current, jdifr , in Eqn. 1. Cur- 
vature dependent sputtering and mass redistribution are 
effectively smoothing at these conditions and can only 
reduce the instability in the ES current. In Fig. |4^b) 
snapshots of the numerical integration of the continuum 
equation and their corresponding 2D FFTs and 2D angle 
distribution are shown. After an integration time of 400 s 
an isotropic pattern forms exhibiting a characteristic pe- 
riodicity without facets. Facet formation starts around 
800 s and at 2000 s facets are already fully developed. 
The 2D angle distribution reveals distinct peaks at 11° 
in diagonal direction of the simulation grid. The 2D FFT 
shows a circular region already with a slight anisotropy 
along the k^- and k^-axes. Finally, at 10000 s very sharp 
peaks appear in the 2D angle distribution at positions ex- 
pected from the zeros of the ES surface current. The 2D 
FFT now also exhibits clear peaks corresponding to the 
4-fold symmetry of the pattern. The comparison with 
the experiments reveals a remarkable agreement. Fur- 
thermore, the proposed continuum equation is able to 
describe the different temperature regimes identified in 
Fig. [I] by choosing the proper coefficients [35] . 

From the experimental fluence series in Fig.[4]we deter- 
mined the temporal evolution of the surface topography. 
In Fig. [5|a) the roughness, w rms , is shown as a func- 
tion of ion fluence. A power law fit to the roughness 
reveals a growth exponent /3 = 0.59, which is close to 
the theoretical value of 0.5 for the so called statistical 
growth limit corresponding to growth by random deposi- 
tion [50] . However, due to the existence of an instability 
larger exponents can be expected [36]. Furthermore, the 
pattern coarsens with ion fluence, i.e. the characteristic 
length increases with irradiation time, i.e. fluence. In 
Fig. [5|b) the characteristic length of the pattern, A, de- 
termined from the height-height correlation function of 
the AFM height images is plotted as a function of ion 
fluence. The characteristic length of the pattern, A, in- 
creases according to a power law with a coarsening expo- 
nent 1/z = 0.14. Such small coarsening exponents have 
also been observed for mound coarsening in homoepitaxy 
[TO] , Theoretically a coarsening exponent of 1/4 is pre- 
dicted for an infinite ES barrier [24]. Smaller exponents 
are expected for moderate barriers in step edge diffusion 
[24] . The numerical integration of the continuum equa- 
tion gives a growth exponent of f3 = 0.45 and a coarsen- 
ing exponent of 1/z = 0.20, in fair agreement with the 
experiments. 

Pattern formation by reverse epitaxy is a universal 
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FIG. 5. (a) Roughness and (b) characteristic length of sur- 
face patterns on Ge(100) as a function of ion fluence. The 
characteristic length was extracted from height-height corre- 
lation functions along the slow and fast scan direction of the 
AFM images, respectively. The lines represent power law fits 
to the data. 



mechanism and can be achieved on many different crys- 
talline materials. We identified the temperature window 
where patterns of crystalline structures are formed, i.e. 
the irradiation temperature is (i) above the dynamic re- 
crystallization temperature of the material and (ii) low 
enough to establish an active ES barrier. At these con- 
ditions an excess of vacancies is created which are par- 
tially reflected at terrace steps inducing an effective up- 
hill mass current. This instability leads to the formation 
of periodic patterns of inverse pyramids oriented along 
the crystalline directions of the surface. The faceting 
of the pyramids results from anisotropic surface currents 
due to a kink ES barrier. We could thus demonstrate 
that ion irradiation can induce patterns of facetted crys- 
talline structures by a reverse epitaxy process in anal- 
ogy to epitaxial growth. The formation of such ordered 
nanostructured surfaces in reverse epitaxy induced by ion 
irradiation is considerably easier and less delicate than 
mound formation in molecular beam epitaxy, where per- 
fect surface preparation and ultrahigh vacuum conditions 
are crucial. Therefore, this technique could establish as 
a complementary epitaxial method for the fabrication of 
high-quality crystalline semiconductor nanostructures. 
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